Occurrence of exponential relaxation in closed quantum systems 
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We investigate the occurrence of exponential relaxation in a certain class of closed, finite systems 
on the basis of a time-convolutionless (TCL) projection operator expansion for a specific class of 
initial states with vanishing inhomogeneity. It turns out that exponential behavior is to be expected 
only if the leading order predicts the standard separation of timescales and if, furthermore, all higher 
orders remain negligible for the full relaxation time. The latter, however, is shown to depend not only 
on the perturbation (interaction) strength, but also crucially on the structure of the perturbation 
matrix. It is shown that perturbations yielding exponential relaxation have to fulfill certain criteria, 
one of which relates to the so-called "Van Hove structure". All our results are verified by the 
numerical integration of the full time-dependent Schrodinger equation. 

PACS numbers: 05.30.-d, 03.65.Yz, 05.70.Ln 



I. INTRODUCTION 

A substantial part of linear non-equilibrium thermody- 
namics essentially relies on a description by means of rate 
equations, often in the form of master equations [![. The 
crucial quantities, such as the probability to find the sys- 
tem in some state i, j or the amount of particles, energy, 
etc. at points i,j in some space, are routinely believed 
to follow equations like 

3 j 

with time-independent transition rates from i to j, R(i — > 
j). Pertinent examples are the decay of excitations in 
atoms, nuclear decay, etc. But also diffusive transport 
phenomena belong to that class, since the diffusion equa- 
tion can also be formulated to take the above form (ran- 
dom walk dynamics). Another implementation of that 
scheme is the (linear) Boltzmann equation [3, [3] where 
particle scattering is taken into account by means of tran- 
sition rates, and many more could be named. 
However, regardless of the incontestable success of such 
descriptions, the strict derivation of rate equations from 
underlaying principles often remains a problem. Typ- 
ically, the descriptiveness by means of rate equations is 
taken for granted. Since those rate equations yield an ex- 
ponential decay towards equilibrium, the basic question 
may be formulated as: How can an exponential decay of 
some observable be derived from the Schrodinger equa- 
tion? 

On the basis of quantum mechanics the most popu- 
lar approach to this question is probably Fermi's Golden 
Rule Despite the undisputed descriptive success of 
this scheme, it is simply derived from first order pertur- 
bation theory, e.g., its validity generally breaks down on 
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a timescale much shorter than the resulting relaxation 
time. Therefore it can hardly describe a complete decay 
into equilibrium. One of the few concrete, concise deriva- 
tions of exponential decay is the Weisskopf-Wigner the- 
ory for the relaxation of excitations in an atom due to the 
coupling of the atom to a zero-temperature, broad-band 
electromagnetic field Q. However, this theory is hardly 
generalizable, since it only applies if just one state is cou- 
pled to a multitude of others, rather than many states 
coupled to many others, as is typically the case. 
A more abstract, rather fundamental approach has been 
suggested by Van Hove [f| It is based on (infinite) 
quantum systems having continuous state densities and 
interactions which are described by smooth functions 
rather than discrete matrices. However, a lot of the find- 
ings for discrete systems in the paper at hand are quite 
parallel to Van Hove's, as will be pointed out below. 
Other approaches are based on projection operator tech- 
niques, in particular the well-known Nakajima-Zwanzig 
(NZ) method. This method is commonly used in the 
context of open quantum systems, i.e., systems that al- 
low for a partition according to a considered system (or 
simply "system" ) and an environment [l|, 0] ■ For a spe- 
cific choice of the initial condition, as pointed out below, 
the projection onto the system's degrees of freedom even- 
tually leads to an autonomous master equation describ- 
ing the dynamics of the system, based on a systematic 
perturbation expansion. But in general, due to the com- 
plexity of higher orders, only the leading order is taken 
into account. In the paper at hand we will demonstrate 
that this truncation may produce wrong results even and 
especially for the case of fast decaying correlation func- 
tions and arbitrarily weak interactions. 
A further approach to this topic is based on the descrip- 
tion of quasi-particle dynamics in many-particle systems 
by the use of Green's functions Q . These considerations 
indicate the validity of a Boltzmann equation. 
In the present paper we will employ another projection 
operator technique-the so-called time-convolutionless 
(TCL) method % M, HO, E3, 0, H- In the follow- 
ing we will follow the TCL-method as detailed in [l4l |. 
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In Sec. [n] we introduce our rather abstract Hamilto- 
nian for a "closed quantum" system (consisting of an 
unperturbed part and a perturbation) and define an also 
rather abstract observable, the dynamics of which we are 
going to investigate. In Sec. IIVI we demonstrate how 
the TCL technique can be used to compute the above 
dynamics of the variable. (This is somewhat reminis- 
cent of projection techniques using "correlated projec- 
tors" [la I16J.) We tune our models such that a lead- 
ing order truncation predicts exponential decay. For a 
"random interaction" this prediction turns out to be cor- 
rect, as is verified by the numerically exact solution of 
the full time-dependent Schrodinger equation. In the fol- 
lowing Sec. W\ non-random ("structured") perturbation 
matrices are discussed in more detail. While a leading 
order truncation still predicts exponential relaxation, it 
is demonstrated that this prediction may fail even for ar- 
bitrarily large models and arbitrarily small interactions. 
This breakdown stems from the fact that higher order 
contributions are not negligible if the interaction matrix 
violates certain criteria. Before we will close with a sum- 
mary and conclusion in Sec. lVIl these criteria will be also 
related to those conditions which Van Hove postulated in 
order to explain the occurrence of exponential relaxation. 



II. MODELS, OBSERVABLES AND 
INTERPRETATION OF DYNAMICS 



In the present paper we will analyze quantum mod- 
els which are much simpler than most of the examples 
mentioned in the introduction. They are defined on a 
very general, rather formal level and are not meant to 
describe any specific, realistic quantum system in great 
detail. The Hamiltonian is taken to consist of a local part 
Hq and an interaction part V such that H = Hq + V . In 
particular, V is assumed to take the special form of an 
"off-diagonal block structure" in the eigenbasis of Hq, 
that is, the matrix representation of H may be written 
as 



H 



■Se 



n — 1 



(2) 



or, equivalently to the above notation, H = Hq + V may 
also be written as 



Ho = Se i*x*i + E -^r 5e \m\ 

n — 1 n — 1 



i=0 
n-1 



3=0 



V = [ v * N)01+H.c. ] , 



(3) 



where \i), \ j) form the basis in which (|2|) is represented. 
Obviously, the complete Hilbert space is divided into two 
subspaces, where i runs through the states of the first 
and j through the states of the second subspace, respec- 
tively. Obviously, Hq may correspondingly be separated 
into two parts which we only specify very roughly at 
this point by two parameters: There are two identical 
"bands" with width Se and n equidistant energy levels 
each. 

The average strength of the interaction V is measured by 



-. n—1 

a 2 = ^EM 



ir 



(4) 



i,j=0 



In our first example in Sec. IIVI we take the matrix ele- 
ments Vij in the off-diagonal blocks to be Gaussian, com- 
plex, random numbers. For the other examples V will be 
specified below. In all cases the matrix elements of V in 
the diagonal blocks are all zero, just to keep the picture 
as simple as possible. 

We will investigate the (relaxation) dynamics of an ab- 
stract observable a, represented by an operator A, which 
is chosen in such a way that 



[A, H ]=0, Tr{A} = , Tr{A 2 } = 1 . 



(5) 



The first of these properties states that A is diagonal 
in the eigenbasis of Hq, while the remaining two prop- 
erties do not mean crucial restrictions on A. While 
all of the following will be correct for any A featuring 
the above properties, we mainly concentrate in our ex- 
amples on "binary" operators, i.e., operators featuring 
only two different eigenvalues, namely, + l/v2 n in one 
subspace and — l/\2n in the other. This means that 
a = Tr{Ap} = + l/^/2n indicates that the system en- 
tirely occupies one subspace and a = — l/y/2 n indicates 
that it entirely occupies the other subspace. (Here, p is 
the density matrix for the state of the system.) If and 
only if a(t) is found to relax exponentially to zero, the 
system allows for a merely statistical interpretation en- 
tirely beyond quantum physics: it is then in accord with 
a system featuring two distinguishable states in between 
it can "hop" with a given transition rate, the latter being 
equal for both directions, a then represents the difference 
between the probabilities of finding it in one or the other 
state, respectively. 

In an abstract way the above model may represent many 
physical situations. It may be viewed as a simplified 
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model for the exchange of an excitation between, e.g., 
two weakly coupled atoms, molecules, quantum dots, etc. 
A then represents the probability to find atom 1 excited, 
subtracted by the probability to find atom 2 excited, V 
represents the coupling in this scenario. Or it may model 
the momentum dynamics of a particle bound to one di- 
mension which possibly changes its direction (forward- 
backward) due to some scattering. In a many-particle 
system the current operator could be identified with A 
and V may stand for a particle-particle interaction. This 
way the dynamics of the current autocorrelation func- 
tion could be investigated based on the framework be- 
low. More detailed information about such models can 
be found in 0, 0, H HI M, HI • 

III. TCL SCHEME AND CHOICE OF THE 
PROJECTION OPERATOR 

In this section we give a short overview of the time- 
convolutionless (TCL) projection operator technique [H, 
[T3j | . Furthermore, we introduce the pertinent equations 
which are applied to models with various interactions in 
Sec. IIVI and Sec. [V] A detailed derivation of these equa- 
tions is beyond the scope of this paper and can be found 

in qui. 

The TCL method is a projection operator technique such 
as the well-known Nakajima-Zwanzig technique [23L |24|. 
Both are applied in order to describe the reduced dynam- 
ics of a quantum system with a Hamiltonian of the type 
H = Hq + V . Generally, the full dynamics of the system 
are given by the Liouville-von Neumann equation, 

^p(t) = - l [V(t),p(t)}=C(t)p(t). (6) 

(Now and in the following all equations are denoted in 
the interaction picture.) In order to describe the reduced 
dynamics of the system, one has to construct a suitable 
projection operator V which projects onto the relevant 
part of the density matrix pit). V has to satisfy the 
property V 2 p(t) = V p(t). Recall that in our case the 
relevant variable is chosen as the expectation value a(t) 
of the binary operator A. For initial states p(Q) with 

Vp(0)=p(0) (7) 

the TCL method yields a closed time-local equation for 
the dynamics of V pit), 

j t Vp{t) = Kit)V P {t) (8) 

with 

oo 

fc(t) = ££i(t). (9) 

The TCL technique avoids the usually troublesome time 
convolution which appears, e.g., in the context of the 



Nakajima-Zwanzig technique. Eq. ((8j) and ([9]) represent 
a formally exact perturbative expansion. 
A brief comment on initial conditions should be made 
here. If ([7]) is not fulfilled, of course an additional inho- 
mogeneity appears on the r.h.s. of ([8]). This may change 
the solutions of ((5]) drastically, c.f. [25[ and references 
therein. However, for the model to be addressed below, 
there is substantial numerical evidence that, for a large 
set of initial states that do not fulfill ||7J), the dynamics 
are nevertheless reasonably well described by ([5]) (with- 
out inhomogeneity) [H HE HE US H [H[ . Having men- 
tioned this issue we consider in the following exclusively 
initial states in accord with 0. 

For many models the odd cumulants of the expansion 
(|9"|) vanish: K 2 %+i{t) = 0. This will turn out to apply to 
our model as well. The lowest non-vanishing order scales 
quadratically with A and reads 

JC 2 (t)= [ dt 1 VC(t)C(t 1 )V . (10) 
Jo 

For the fourth order term one finds 
pt pt\ pt 2 
JC 4 (t) = dh dt 2 / dt 3 
Jo Jo Jo 

V Cit) dh) Cih) C{t 3 )V 

- VC{t)C{h)VC(t 2 )C{t 3 )V 

- vc{t)c(t 2 )rcih)c{t 3 )r 

- VC{t)C{t 3 )VL{ti)L{t 2 )V . ill) 

Note that the TCL approach is commonly used in the 
context of open quantum systems [l], 0, HE H3| . The 
TCL method is, however, also applicable to our closed 
quantum system. 

To those ends, we define the projection operator V by 

v p{t) = h 1 + A Tr{ A p(t)} = h l+A a[t) ■ (12) 

As already mentioned above, V is constructed to project 
onto the time-dependent expectation value a(t) of the 
binary operator A, in the Schrodinger picture. But since 
A commutes with Hq, this expectation value is identical 
in the interaction and the Schrodinger picture. The full 
dynamics [Hilbcrt space: dimension 2n, Liouvillc space 
of density matrices: dimension (2n) 2 ] is broken down to 
the time evolution of the single variable ait), all other 
information is neglected. As a suitable initial condition 
we can then choose p(0) = (l/2n) 1 + (1/V2n) A which 
implies a(0) = l/V2n. Inserting Eq. (12]) into Eq. © 
yields the closed equation 

oo 

a{t) = Y J K i {t)<t) (13) 

i=l 

with Ki{t) = Tk{AKi{t) A}. Due to Eq. (jT2j>, the second 
order term reads 

K 2 (t) = - f dt'C(t') , (14) 
Jo 



4 



where the two-point correlation function C(t') is given 
by 

C{t') = Tt{i[V{t),A] i[V(ti), A]\ , if = t-ti . (15) 

A rather lengthy but straightforward calculation yields 
for the fourth order 





- 




Jo 


h = 


H 


h = 


- C(t 


h = 


- C(t 


h = 


- C(t 



dti 



ti 



dt 9 



dt 3 h + h + h + h 



(16) 



IV. SECOND ORDER TCL AND COMPLETELY 
RANDOM INTERACTION 

In this section we apply the equations in second order 
TCL to a model with the completely random interaction 
introduced in Sec. HU The function C(t') in Eq. (fTB")) is 
identical to the autocorrelation function of the interac- 
tion, since it can also be written as 



C(t>) 



^ \ V lj\ 2 COS[Wy {t-tx)] 



(17) 



i,j=0 



with frequencies w^j = (i—j)/(n—l) Se corresponding to 
Hq. Here, just like in many other examples, C(t') de- 
cays within the correlation time tq which is of the order 
of tc ~ 4ir/Se for our model. Afterwards the integral 
K2(t) becomes approximately time-independent and as- 
sumes a constant value R until the "Heisenberg time" 
T = 2irn/Se is reached. This behavior can be inferred 
from integrating (|17p and exploiting the properties of the 
sinc-function. From this analysis also R may be found 
with an accuracy determined by the law of large num- 
bers. Thus the second order approximation of Eq. OX 
eventually results in 



a(t) = —R a(t) , R . 



47rn A 2 
Se 



(18) 



We hence obtain a rate equation featuring the form of 
Eq. and thus exponential dynamics for a(t). The 
solutions for a(t) decay exponentially with a relaxation 
time Tr = 1/R. However, this result is only valid within 
the boundaries tq <C tr < T, because K 2 {t) can only 
be considered as time-independent up to the Heisenberg 
time. Recall to this end that our model features equidis- 
tant energies such that C(t') is strictly periodic with T. 
These two boundaries also result in two necessary criteria 
for the system parameters which have to be fulfilled in 
order to produce the occurrence of exponential dynamics, 



16vr 2 n A 2 
5? 



< 1 



in 2 n 2 A 2 



Se 2 



> 1 



(19) 



Remarkably, the whole derivation of the rate equation 
using second order TCL does not depend on the details 
of the interaction, i.e., the individual absolute values of 
the single matrix elements as well as their relative phases 
are not relevant. We should already mention here that 
the "structure" , which we are going to introduce into the 
interaction in the following section, only concerns those 
details, hence the second order contribution Ki will be 
the same in all our following examples. 
In Fig. [T] the numerical solution of the Schrbdinger equa- 
tion is shown for the above repeatedly mentioned ran- 
dom interaction and compared with the TCL prediction. 
All parameters (the width of the Gaussian distribution 
according to which the matrix elements of V are gen- 
erated, the bandwith, etc.) are adjusted such that the 
criteria are well satisfied. This solution is obtained 
by exact diagonalization. And in fact, we find a very 
good agreement with the theoretical prediction of second 
order TCL. 




4000 



time t 



FIG. 1: Time evolution of the expectation value a(t) for an 
interaction with completely random Vij . The numerical result 
(crosses) indicates exponential behavior and is in very good 
agreement with the theoretical prediction (|18[1 of second order 
TCL (continuous curve). The system parameters n — 1000, 
Se = 0.5, A = 2.5 • 10~ 4 fulfill the conditions $P3$. 



V. FOURTH ORDER TCL AND NON-RANDOM 
INTERACTIONS 

In this section we will be concerned with the structure 
of the interaction matrix and, especially, its influence on 
the time evolution of the expectation value a(t). It will 
be demonstrated that the theoretical prediction (fT5)) of 
second order TCL fails to describe the numerically exact 
solution of the Schrodinger equation correctly for certain 
"interaction types" , even and especially if the conditions 
(UHJ) are fulfilled, i.e., the "strength" is "adequate". We 
will outline that this failure stems from the fact that the 
fourth order contribution of the TCL expansion is not 
negligible on the relaxation timescale which is obtained 
from second order TCL. However, the exact evaluation 
of Ki(t) turns out to be almost impossible, analytically 
and numerically. Instead we will present feasible estima- 
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tions of Ki(t) I K<z(t) based on suitable approximations of 
K 4 (t) called S{t) (see (|23|38)). Whenever 



s(t) 

K2{t) 



< 1 



(20) 



is violated the influence of higher order terms is not neg- 
ligible. If this is the case for times t of the order of or 
shorter than tr, no exponential relaxation will result. 



A. Uniform Interactions and Van Hove structure 

Let us start with an example. Fig. [2] shows the time 
evolution of the expectation value a(t) for an interac- 
tion with Vij — A. This type of interaction is, of course, 
highly non-random, since all matrix elements have the 
same absolute value and phase. The second order ap- 
proximation obviously yields a wrong description for this 
interaction structure, that is, the dynamics are not ex- 
ponential, although both of the conditions p^|) are well 
fulfilled. It should be remarked again that the observed 
non-exponential behavior definitely is a structural issue. 
For instance, it can not be "repaired" by simply decreas- 
ing the overall interaction strength, because this decrease 
would eventually lead to the violation of the criteria (fT9)) . 



e 




25000 



time t 



FIG. 



2: Time evolution of the expectation value a(t) for an 
interaction with Vij = A. The theoretical prediction (|18p of 
second order TCL (continuous curve) fails to describe the 
numerical solution (crosses) correctly, although the system 
parameters n = 1000, 5e = 0.5, A = 2.5 • 10" 4 still fulfill the 
conditions (1191). V violates the Van Hove structure. 



To analyze this model we now develop our first estimate 
S(t) for K 4 {t) which concerns the timescale t >~ tc- We 
start from from Eq. (|16[l . where we abbreviate the triple 
time integration by a single " f " . One may hence write 
K i{t) = I h + h + h + h- Fig. [3] shows a sketch for the 
integration volume of K 4 (t) in the 3-dimensional space 
which is spanned by t\ 1 t^. The integration does not 
run over the whole cube with the edge length t, but only 
over the region where t% < ti < t\ holds. 
C(t—tx) is the autocorrelation function of the interaction 
which has already been mentioned in (|15p . Recall that 




FIG. 3: Sketch for the integration volume of K&(i) for fixed 
t. The cube with the edge length t is drawn with thin lines, 
the actual integration volume is marked with thick lines. The 
dashed lines represent the changed integration volume which 
is used in the approximations. 



C(t') is only different from zero around t' = in a small 
interval of the width tq- Thus, the integrands I2, I3, I a 
are only different from zero in a small volume around the 
region where both of the arguments are equal for each of 
the two multiplied correlation functions. 
First of all let us focus on I3 as well as I4. The integrand 
J 3 contributes to K 4 {t) for t = ti and t\ = £3, while the 
integrand I4 contributes to K 4 {t) for t = £3 and t\ =ti-, 
respectively. The sketch in Fig. [3] displays that both of 
these regions overlap only in the vicinity of one single 
point with the integration volume of K 4 (t), namely, at 
the point where all arguments are equal to t. Especially, 
this overlap does not increase with t. Therefore the triple 
time integration is estimated by f h ~ J I 4 w C(0) 2 Tq. 
Using the estimate R w C(0)tc, we eventually obtain 
for the ratio between the contributions from J3 , 14 to the 
fourth order and the second order Kz(t > re) for times 
t > T C 



III 
R 



[h 
R 



c(oy 



TR 



(21) 



Recall that the derivation of exponential behavior within 
second order TCL has required tq <C tr or, equivalently, 
ct« 1 such that the contributions to K 4 (t) which arise 
from I3 and 1 4 are negligible, at least in comparison with 
R. 

Analogous conclusions cannot be made for the term I2, 
because its overlap with the integration volume is larger 
and grows with t. We have to find another estimation 
for the contributions of I2 as well as I\, of course. Our 
estimation is based on the fact that neither I2 nor I± 
can decay on a shorter timescale than tq in any possible 
direction of the (ti, t 2 , ts)-space. This fact is obviously 
correct for 12- But what about Ji? Since the term 1\ 
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consists of summands which have the typical form 

v ab v bc v cd v da e~ l 4 e" 4 u e~ l * 2 e" 1 * 3 , (22) 

only those frequencies and, especially, those largest fre- 
quencies in Ii which have already appeared in C(t') con- 
tribute significantly to I±. Consequently, I\ can never 
decay faster than C(t') in any possible direction of the 
(ti,t2,^3)-space. In the (possibly unrealistic) "best case" 
I\ + 1 2 decays within t c around the point ii/ 2 (i, t, t, t) — 
^1/2(0,0,0,0) = /i/2(0). We can therefore estimate the 
value of q(jc) by 



|/i(Q)-c(Q) 2 ]^ 



R 



C(0)t c 



h(0) 



(23) 

[3 is a lower bound for the ratio between the fourth and 
the second order of TCL for times t > tc- If (3 ~ 1 or 
even larger, then K±{t > Tc) dominates Ki(t > Tc), that 
is, exponential behavior in terms of the second order 
prediction cannot occur. But j3 -C 1, however, does not 
allow for a strict conclusion, since a slower decay of I\ or 
I2, as the case may be, raises their contribution to Ki{t). 
Nevertheless, the condition (3 <C 1 is an additional 
criterion for the occurrence of exponential decay which 
involves the structure of V. 

In the following we will discuss why and to what ex- 
tend (3 and, especially, the ratio /i(0)/G(0) 2 is related to 
the conditions which have been postulated by Van Hove 
for the interaction V in order to explain the onset of ex- 
ponential relaxation, see To this end, let us define 
a hermitian operator G by 

G=[V,[V,A}}. (24) 

A straightforward calculation yields 

/ 1 (0) = Tr{G 2 } = ^|G y | 2 , G(0) = Tr{AG}, (25) 

where G^ represents the matrix elements of G in the 
eigenbasis of H . Furthermore, let us also introduce the 
superoperator T> which is given by 



VM = J2\i)Mi, 



(26) 



and projects any operator M onto its diagonal elements 
in the eigenbasis of Hq. Then the expression 



(A,G) = Tr {AVG} 



(27) 



defines an inner product between the operators A and G, 
because (A, G) = (G,A)*[= (G, A)] holds and (A, A) = 
1, as well as (G, G) = JT^ G 2 i: are both positive, real num- 
bers. The Schwartz inequality (A, Gf < (A, A) 2 (G, G) 2 
can consequently be formulated. By the use of (A, G) = 
G(0) we eventually obtain 



c(o) 3 <£G&<£l G tfl 2 = -M°) 



(28) 



i.e., ii(0)/G(0) 2 > 1. G(0) 2 is at most as large as the sum 
of the squared diagonal elements of G, according to the 
above equation. Therefore ii(0)/G(0) 2 « 1 and hence 
sufficiently small (3 can only be realized if the diagonal 
elements of G and thus the diagonal elements of V 2 are 
as large as possible in comparison with the remaining 
non-diagonal elements of V 2 (G). In principle, this is 
essentially what Van Hove proclaimed [g, Q . 
In this sense, we define the "Van Hove structure" in the 
context of finite quantum systems: The interaction V is 
said to feature Van Hove structure if 



m 

G(0) 2 



a < 1 



(29) 



while all conditions of second order TCL are simultane- 
ously kept, of course. The latter refers to the validity 
of Eq. (p~9|) . The comparison with (|23|) shows that the 
Van Hove structure implies j3 <C 1 and hence the relax- 
ation may possibly be exponential, as described by the 
second order. Since the evaluation of (3' is much more 
efficient than the complete computation of fourth order 
TCL (there is no time dependence left, e.g., ii(0) only 
depends on t = t\ = f 2 = £3 = 0), the Van Hove struc- 
ture eventually is an assessable criterion for the possible 
occurrence of exponential decay. It is a criterion in the 
sense that only if (|29|) is satisfied, a use of the second 
order approximation is justified for any time longer than 
the correlation time, i.e., t > tq- 

Let us now apply these results to the already introduced 
models with random and non-random (vij = A) interac- 
tions, respectively. The only term which varies for the 
different models is /i(0), since the terms G(0) 2 w n 2 X 4 
and a w 16ir 2 nA 2 /8e 2 (again with an accuracy set by the 
law of large numbers for the random interaction) are the 
same for random and non-random interactions. For the 
random interaction a straightforward calculation leads to 



h(0) =32n 2 A 4 , /3' = 2a<l 



(30) 



such that the random interaction indeed features Van 
Hove structure. This agrees with the numerical results 
in Fig. [1] which yielded exponential relaxation. In the 
case 'Vij = A, however, we finally obtain 



7i(0) = 16n 3 A 4 , (3' = 



16ir 2 n 2 A 2 
5? 



(31) 



where ft > 1, according to Eq. (fl"9]) . The absence of the 
Van Hove structure already suffices to explain the break- 
down of exponential behavior in Fig. [2l 
One may nevertheless be inclined to argue that the Van 
Hove structure is not the crucial difference between those 
two cases but simply the randomness of the matrix el- 
ements (which possibly induces quantum chaos). We 
therefore present a counter-example which immediately 
disproves such an argument. The example is slightly dif- 
ferent from the others, since the complete system is not 
partitioned into equally large subspaces. n\ and 712 de- 
fine the number of levels of the respective subspaces. One 
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subspace consists of only one state (ni = 1). Thus, in 
the matrix V there is only a single column with non- 
zero elements and a single row, respectively. Although 
these non-zero elements are chosen to be all equal (non- 
random) , it can be shown that this V features Van Hove 
structure. Note that such a Hamiltonian occurs, e.g., in 
the context of spin-boson models at zero temperature or 
the scenario addressed by the Weisskopf-Wigner theory, 
see 01. 




4000 
time t 

FIG. 4: Time evolution of the expectation value a(t) for the 
interaction of spin-boson type. The numerical result (crosses) 
indicates exponential behavior and perfectly agrees with the 
theoretical prediction (|18[1 of second order TCL (continuous 
curve). System parameters: n\ = 1 (single level), 722 = 2000 
(many levels), 5e = 0.5, A = 2.5 ■ 10" 4 . 

Fig. 2] shows an almost perfect correspondence between 
the numerical solution of the Schrodinger equation and 
the theoretical prediction (fT8|) which is obtained by the 
use of second order TCL. Here, exponential relaxation is 
found, although V is not randomly chosen. 

B. Sparse Interaction and Localization 

So far, we numerically found exponential decay in ac- 
cord with the second order for all considered models that 
showed the Van Hove structure. There is, however, non- 
exponential behavior for some types of interactions which 
feature the Van Hove property in the sense of (|29|) and 
are in accord with (fll)]) . Recall that those are only nec- 
essary but not sufficient conditions for the occurrence of 
exponential decay. 

An example for such a situation is a model with a random 
but, say, "sparsely populated" interaction. This model is 
almost identical to the model with the completely ran- 
dom interaction. The only difference is that only 1/10 of 
the matrix elements are Gaussian distributed numbers, 
all others are zero. The non-zero numbers are randomly 
placed. Apparently, this type of interaction fulfills the 
Van Hove structure, since the completely random inter- 
action already does. 

Fig. [5] displays the numerical solution of the Schrodinger 
equation and the theoretical prediction (fTS)) of second 
order TCL. At the beginning there is a good agreement 








4000 
time t 

FIG. 5: Time evolution of the expectation value a(t) for an 
interaction with random but "sparsely populated" Vij. The 
theoretical prediction (|18|l of second order TCL (continuous 
curve) deviates from the numerical solution (crosses), even 
though the Van Hove structure as well as the conditions (|19p 
are fulfilled. System parameters: n = 1000, Se — 0.5, A = 
2.5 ■ 10~ 4 . 



but then the numerical solution starts to deviate from a 
purely exponential decay and finally sticks at a clearly 
positive value. The latter non-zero value may be a hint 
towards localization effects which also appea r, e.g., in 
the context of the Anderson model [3(1 l3lll32l. l33j . And 
in fact, the sparsely populated interaction takes a form 
which is very similar to the Hamiltonian of the, e.g., 3- 
dimensional Anderson model in the chaotic regime. 
Apparently, we have to extend the analysis of the fourth 
order: There is no exponential behavior by the means of 
a complete exponential decay, although V fulfills the Van 
Hove property. Recall that the Van Hove criterion has 
been derived from the consideration of times t tq 
and thus t = t\ = ti = £3 <« re- Hence, we have to 
reconsider the full time dependence of the fourth order 
to produce a feasible estimate for the timescale t s» tr. 
To this end, the integrand I\ is expressed by 

h =Tr{G{h,t)G(t 2 ,t 3 )} , (32) 

where the hermitian operator G(t±,t) is again given by 

G(t,ti) = [V(t),[V(ti),A]]. (33) 

If li(0) « C'(O) 2 , the diagonal terms dominate at t = 
ti = t% = £3. Based on this fact, we carefully assume 
that 1\ is dominated by these terms for other times as 
well. Of course, this assumption neglects the larger part 
of all terms but leads, as will be demonstrated below, to 
a criterion which may be evaluated with limited compu- 
tational power. (For our simple example its validity can 
also be counterchecked by direct numerics.) However, 
following this assumption, 1\ can be approximated by 

JiM^GwCt-tOGiiCta-ts), (34) 

i 
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where Gu(t — t\) are the diagonal matrix elements of 
G(t,ti) in the eigenbasis of H , namely, 

Guit-h) =2Y^{A ii -A jj ) \V %3 \ 2 cos^f-fi)] . (35) 



Furthermore, the correlation function C'(t — t\) can, by 
the use of this notation, also be written as 



C(t-t 1 ) = Y,A ii G ii (t-t 1 ) 



(36) 



such that 12, the remaining fourth order integrand, can 
be expressed as well by 



i,3 



G u {t 



t{) Ajj Gjj(t 2 - h) 



(37) 



In order to estimate with reasonable computational effort 
how Ki(t) compares with Ki(t) another approximation is 
necessary. Obviously, the expressions for 1\ , I2 are invari- 
ant along lines described by t\ — const., £2 = £3. Thus, as 
an approximation, we shift the integration volume from 
the original region, indicated with solid lines in Fig. [3J to 
a new region, indicated with dashed lines in Fig. [3] Obvi- 
ously, this is a rather rough estimate but it will turn out 
to be good enough for our purposes. Now the coordinate 
transformation x = t—ti, y = £2 — ^3, z = t — t2 decouples 
the integrations within the new integration volume such 
that we eventually find for S(t) « K±(t) (if V features 
Van Hove structure) 



S{t)=t[£T l {t?-K2{tf 



(38) 



with the time integral T % (t) = f* dt' G u (t'). Now 
may eventually be checked with very low computational 
power, based on S(t) from (|38p . This adds to fT9")) and 
further manageable criterion for exponential 
relaxation. Fig. [5] shows q(t) [based on (f3"5| ] for the fol- 
lowing interaction types: the completely random interac- 
tion, the interaction of spin-boson type, and the random 
but sparsely populated interaction. Fig. [5] apparently 
demonstrates that this approximation is able to explain 
the breakdown of exponential behavior in the case of a 
random, sparsely populated interaction: The fourth or- 
der becomes roughly as large as the second order at a 
time which agrees with the deviation between the second 
order theory and the numerical results in Fig. 03 In both 
other cases q(t) remains sufficiently small, at least until 
the relaxation time is reached. 

Obviously, regardless of the interaction type, K±(t) will 
eventually dominate K%{() for large enough times (c.f. 
34]). This, however, does not necessarily spoil the expo- 
nential decay: If a(t) has already decayed almost com- 
pletely into equilibrium, even a significant change of the 
rate K(t) will not change the overall picture of an expo- 
nential decay (, as long as K{t) remains negative). The 
influence of a large K±(t) will only be visible if it occurs, 



<x x 



2000 



1000 
time t 



FIG. 6: Time evolution of the value q(t) based on (|38|) . q(t) 
is numerically calculated for a completely random interaction 
(crosses), the interaction of spin-boson type (squares), and 
the random but sparsely populated interaction (circles) . Note 
that Kz(t) is a constant which is also identical for all three 
interactions. The system parameters are chosen according to 
Fig. [1] Fig. |4] and Fig. [5] respectively. 



while a(t) is still far from equilibrium, i.e., at times of 
the order of tr. If one now computes the ratio q(t) for 
the time t = tr, one finds 



q(tr) 



S,r,(r fl ) 2 

R 2 



1 



(39) 



where one has to take (|3"5]l and K 2 (t r ) = R = 1/tr 
into account. This form has the advantage of being com- 
pletely independent of the overall interaction strength A. 
One can hence compute q(Tu), taking tr as a free vari- 
able. The region in which q(rn) <« 1 then represents 
the range of different tr for which exponential decay is 
possible and to be expected. The different "possible" 
tr can then be implemented by tuning A appropriately. 
Often q(r r) is found to increase monotonously, essen- 
tially like in Fig. [Sj Thus a good number to characterize 
a class of models with different relaxation times (inter- 
action strengths) would be r max as the largest time for 
which q{T max ) <« 1 holds true. This then indicates the 
largest timescale on which exponential relaxation can still 
be expected. We should note here that we intend to use 
this measure, T max , to investigate transport behavior in 
models of the Anderson-type in a forthcoming paper. 



VI. SUMMARY AND OUTLOOK 

We investigated the dynamics of some expectation val- 
ues for a certain class of closed, finite quantum systems 
by means of the TCL projection operator method. This 
technique yields a perturbation expansion for those dy- 
namics. Taking only the second (leading) order into ac- 
count, we find that the evolution of these expectation 
values may be described by a rate equation, i.e., they 
relax exponentially if certain criteria are fulfilled. Those 
criteria, however, only depend on "rough" parameters 
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like overall interaction strength, bandwidth and density 
of states but not on, e.g., the phases of the interaction 
matrix elements. An adequately computed numerical 
solution of the Schrodinger equation is in accord with 
this leading order result for random interaction matri- 
ces. However, numerics also show that this accordance 
breaks down if one considers non-random interactions, 
even if the above rough criteria are met. This, of course, 
indicates that higher orders are not negligible, depending 
on the structure, not only on the strength of the interac- 
tion. Subsequently, we established a numerically simple 
estimate for the absolute value of the fourth, i.e., the 
next higher order, in comparison to the second, for short 
times. From this approach it can be inferred that the 
fourth order remains negligible at small times if the in- 
teraction features a certain structure which we define as 
Van Hove structure according to [5]. However, numerics 
indicate that for certain interaction structures the fourth 
order may become non-negligible at larger times, thus 
spoiling the exponential relaxation, even if the interac- 
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